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Abstract 

We have searched for the presence of simple P and S-bearing molecules in 
Titan's atmosphere, by looking for the characteristic signatures of phosphine 
and hydrogen sulfide in infrared spectra obtained by Cassini CIRS. As a 
result we have placed the first upper limits on the stratospheric abundances, 
which are 1 ppb (PH3) and 330 ppb (H2S), at the 2-cr significance level. 
Keywords: Titan, Titan, atmosphere, Abundances, atmospheres. 
Atmospheres, Composition, Saturn, satellites 

1. Introduction 

To date, no molecules bearing the light elements phosphorus (P) and 
sulfur (S) have been detected in the atmosphere of Titan by either remote 
sensing or in-situ methods. However, from cosmological considerations both 
P and S must have been present in the icy planetesimals that formed Titan, 
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and also delivered in trace quantities by later impacts. P in the form of PH3 
is found in Jupiter's atmosphere at approximately solar abundance, while on 



Saturn its abundance is around 3x solar (Owen and Encrenaz, 2003). This 



enrichment is in line with core-accretion models, which predict that Saturn 
had a larger ice-to-gas fraction in its formation compared to Jupiter. Sulfur 



is enriched on Jupiter by 2.5 x solar (Niemann et al. , 1998) (and not yet 
detected on Saturn), perhaps due to easy formation of H2S clathrate hydrates 



Owen and Encrenaz (2003). Both P and S should therefore be present in 



Titan's bulk composition at fractions greater than those on Saturn, since the 
ice to gas ratio must have been much higher. 



Fortes et al. ( 2007 ) considered a possible role for sulfur compounds in Ti- 



tan cryovolcanism, suggesting that ammonium sulfate could form a magma 
in Titan's mantle. Plumes of this magma could dissolve methane present in 
crustal clathrates, allowing explosive release to the surface. One prediction 
of this model is that ammonium or other sulfates should be detectable spec- 
troscopically on Titan's surface, as is the case on Europa and Ganymede. 
Such a model would presumably also release trace amounts of sulfur into 



23 Titan's atmosphere, especially in the wake of eruptions. Pasek et al. (2011) 



focused on the role of phosphorus in Titan, arguing that it could be delivered 
both endogenously and exogenously to the surface. Phosphine is efficiently 



26 trapped in clathrate hydrates, and Pasek et al. (2011 ) show that all clathrates 



that trap H2S also trap at least as much PH3. Moreover, they show that PH3 
is more soluble in organic liquids than in water, and therefore any phosphine 
released from melted clathrates could be dissolved in Titan's hydrocarbon 
lakes, and would participate in the hydrocarbon meteorological cycle. 
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31 We have therefore searched the infrared spectrum of Titan for the sig- 

32 natures of the two most hkely carriers of P and S: phosphine (PH3) and 

33 hydrogen sulfide (H2S). Our data is from the Composite Infrared Spectrome- 



34 ter (CIRS, Flasar et aL , 2004), carried onboard the Cassini spacecraft, which 

35 has been making close flybys of Titan since attaining Saturn orbit in mid- 
36 2004. Both gases have signatures in the 8-11 fim (1250-900 cm~^) range, 

37 which is largely free of the hydrocarbon emissions that dominate much of 

38 Titan's infrared spectrum. We first model and subtract out the emissions of 

39 known species including methane (CH4), acetylene (C2H2), ethylene (C2H4) 

40 and deuterated methane (CH3D). We then add H2S and PH3 incrementally 



41 iJ^o our spectral calculation using standard line lists from HITRAN (Rothman 



et al. , 2009), and compare the predictions to the remaining Titan spectrum. 



43 By comparing the model predictions to the instrument noise level at the 

44 1, 2, and 3-cr levels, we have derived the first numerical upper limits on the 

45 prevalence of H2S and PH3 in Titan's stratosphere. We follow a report of our 

46 results with a discussion of the implications for existing models, and finish 

47 with some conclusions about the prospects for future searches. 



48 2. Method 

49 2.1. Instrument and Dataset 

50 The Cassini CIRS instrument is a dual design, comprising mid-infrared 

51 (1400-600 cm^^) and far-infrared (600-10 cm~^) spectrometers, both with a 

52 FWHM (full-width to half-maximum) spectral resolution variable from 0.5- 

53 15.5 cm~^ after Hamming apodization. The mid-IR Michelson spectrometer 

54 employs two 1 x 10 detector arrays: focal plane 3 (FP3, 600-1100 cm~^) 
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55 and focal plane 4 (FP4, 1100-1400 cm~^). Both arrays have square pixels 

56 with fields of view 0.3 mrad across. FP4 has the highest sensitivity (lowest 

57 background noise) and least instrumental interferences, and was consequently 

58 used in this study. 

59 The observations were described in an earlier paper ( [Nixon et al. , 2010). 



60 The spectra were acquired during the SS*'* Titan flyby (T55) of Cassini on 

61 May 22"'^ 2009, in a four-hour period when the spacecraft was at a range 

62 of 116,000-177,000 km. The observation (known as 'MIRLMPAIR'-type) 

63 was targeted at Titan's limb, with the two arrays above and parallel to the 

64 horizon. FP4 was above FP3, with the pixels centered at an altitude of 247 

65 km (0.27 mbar) and spanning approximately 44 km (just under one scale 

66 height) in the vertical direction at the mid-point of the observation. Using 

67 CIRS PAIR mode, all ten detectors simultaneously recorded data, paired into 

68 five receiver channels. A total of 941 pair-mode spectra at 0.5 cm~^ resolution 

69 were selected and then averaged to create a single, high signal-to-noise (S/N) 

70 ratio spectrum. 

71 2.2. Spectral Modeling 

72 The CIRS FP4 limb spectrum was then modeled to remove the signatures 

73 of known gas species. These included CII4 (z/4 band at 1305 cm~^) and 

74 CII3D (z/g band at 1156 cm~^), plus some weaker contributions from the 

75 hydrocarbons C2H4 and C2H2. The modeling closely follows that described 



76 in a recent paper (Nixon et al. , 2012), and is briefiy summarized here. 

77 An initial vertical atmospheric model was created with 100 layers equally 

78 spaced in log pressure from 1.45 bar to 0.05 yubar, based on the temper- 

79 ature profile and major gas abundances (N2, CH4, H2) determined by the 



Huygens probe (Fulchignoni et al. 2005 Niemann et al. 2010). Other gas 



species (C2H4, C2H2) were included with constant vertical abundances in the 
stratosphere, at initial values taken from previous CIRS measurements at 



low latitudes (Coustenis et al. , 2010). Similarly a uniformly mixed (constant 



84 particles/g atmosphere) stratospheric haze was included with optical prop- 



85 erties from Khare et al. (1984). The forward radiative transfer model was 



86 computed using the NEMESIS code (Irwin et al. , 2008) applied to the model 



atmosphere and using the HITRAN gas line atlas (Rothman et al. , 2009), 



ip,nd then convolved with the FP4 detector spatial response shapes (Nixon 



et al.l 2009b) to generate a model spectrum. 



At this point the model was iterated to arrive at an optimum fit to the 
measured spectrum, by adjusting the model temperature profile (at each 
layer) and uniform vertical gas abundances of C2ll2 and C2H4, to minimize 
a cost function similar to a figure of merit. This approach has been 
successfully used to fit the T55 limb spectrum of Titan in previous work 
dNixon et alj |2009a[ |2010[ ). Fig[T] (a) shows the best fit to the data. Note 
that a weak band of propane (z/7 at 1158 cm^^) was not included in the 
model, as a line list for this band has yet to be produced, and shows clearly 
in the data-model residual (Fig [T] (b)). The strong CH4 2/4 band is fitted 
reasonably well, but also shows some residual mismatch, which may be due 
to imperfect knowledge of the underlying haze (continuum) opacity and/or 
wavelength calibration uncertainties. 



102 2.3. Determination of Abundance Upper Limits 

103 To search for PH3 and H2S, we avoided spectral regions that showed non- 
104 random 'noise' residual: especially the CsHg z/7 region (1140-1200 cm~^) and 



105 the CH4 region {u > 1250 cm~^). This permitted us to search for the z/4 

106 band of PH3 around 1120 cm~^, and for a portion of the weak H2S 1^2 band 

107 at 1200-1250 cm~^ (see Fig[T](c) and (d)). The summed hne strengths in 

108 these regions were 1.40 x 10~^^ cm molecule"^ for 1121 hues of PH3 from 

109 1080-1140 cm-^ and 2.2 x 10~20 cm molecule-^ for 130 hues of H2S at 1200- 



110 1250 cm-\ using hne data from HITRAN ( [Rothman et aL||2009[ ). Note that 

111 the PH3 band intensity was almost 100 x that of H2S in the spectral ranges 

112 considered, so we expect 100 x greater sensitivity to PH3 compared to H2S. 

113 These two gas species were added to the best-fit atmospheric model at 

114 fixed trial abundances of 1 ppb with other gases and temperature held con- 
ns stant at the previously retrieved values, and NEMESIS was allowed to re- 

116 trieve a 'best fit' abundance for the new species. In both cases, these retrievals 

117 resulted in no statistically significant improvement to the model x^, hereafter 

118 designated Xo- the minimum x^- We then proceeded to calculate upper lim- 



119 its to the abundances, following the approach of Teanby et al. (2009) and 



Nixon et al. (2010). Starting with very low trial gas abundances in the model. 



121 these were increased incrementally, and at each value the forward spectral 

122 model was calculated, without optimized fitting or inversion, and the data- 

123 model computed. As the trial gas abundances increase, so too does the 



Xo- See Fig. [2| Following [Press et al.| ( [T992| the 1, 2 and 3-ct 



125 upper hmits to the gas abundances occur at Ax^ = 1, 4, 9. 



126 3. Results and Discussion 

127 Table 1 shows our results. The derived 1, 2, 3-cr maximum abundances for 

128 PH3 are 0.3, 1.1, 2.2 ppb respectively, while the corresponding values for H2S 
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129 are 91, 330, 700 ppb. The H2S upper limits are some two orders of magnitude 

130 higher than those of PH3 as expected, due to the much weaker band inten- 

131 sity in the spectral region considered, along with a higher 1-a NESR (Noise 

132 Equivalent Spectral Radiance) level in the corresponding spectral interval 

133 (0.6 versus 0.3 nW cm~^ sr~^/cm~^). 

134 Sulfur compounds have been suggested to play a role in Titan cryovol- 



135 canism (Fortes et al. , 2007), and therefore be present in trace amounts in 

136 Titan's atmosphere. Phosphorus is also easily dissolved in organic solvents, 

137 and as such could be a component of Titan's present-day methane-ethane 



138 lakes, and participate in the hydrocarbon cycle (Pasek et al. , 2011). A sim- 

139 pie calculation tracing the saturation vapor pressure (data from public NIST 

140 and CRC databases) up through Titan's troposphere using the Huygens tem- 

141 perature profile, indicates that abundances of PH3 and H2S that could reach 

142 the stratosphere are 6 and 0.035 ppb respectively. Therefore, our 2-a upper 

143 limit of 1 ppb for PH3 provides some constraint on tropospheric PH3, while 

144 our H2S limit of 330 ppb is four orders of magnitude greater than allowed in 

145 this scenario and is not constraining. 

146 At the present time, the lack of detection of P or S-bearing species any- 

147 where in Titan's atmosphere has caused these elements to be omitted from 

148 present photochemical models, so we have no predictions for comparison. In 

149 future, it could be interesting for such models to allow for some injection 

150 of P and S molecules into the atmosphere through a putative eruption, and 

151 to investigate subsequent chemistry, for example, conversion of sulfates to 

152 sulfides. 
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153 4. Conclusions and Further Work 



154 In this paper we have searched for simple compounds of phosphorus and 

155 sulfur in Titan's infrared spectrum recorded by Cassini CIRS, placing the 

156 first upper limits on the abundance of PH3 and H2S. In the stratosphere 

157 at 247 km we find that no more than 1 ppb of PH3 and 330 ppb of H2S 

158 can be present, at the 1-a level of significance. Some potential exists for 

159 improving on these upper limits using CIRS, e.g. by future measurements 

160 at a lower limb altitude where the atmospheric density is greater. The peak 

161 sensitivity for weak trace gases is often near 10 mbar (~100 km), at least at 

162 low latitudes, so a repeat measurement in the lower stratosphere may yield 

163 more stringent limits. 

164 Searching in the far-IR spectrum may prove helpful in the case of H2S, 

165 which has rotational lines at 50-150 cm~^ that are some 100 x stronger than 

166 the z/2 considered here. However in this region the NESR of CIRS FPl is 

167 more than 10 x higher than that of the FP4 detectors used in this work, and 

168 suffers from systematic noise artifacts. Nevertheless, the sub-mm range has 

169 already proved productive for molecular line searches by instruments such 

170 as Herschel and ground-based sub-mm telescopes, resulting in the detections 



171 of CH3CN dMarten et aLf |2002D and HNC ( [Moreno et sX\ [2011j), and wiU 

172 doubtless reveal further new species in due course. 
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Table 1: Calculated abundance upper limits for PH3 and H2S 







Pressure 




Wavenumber 


l-(7 


Upper Limits (ppb) 


Gas 


Lat. 


(mbar) 


Band 


Range (cm~^) 


NESR* 


1-0" 2-0" 3-0" 


PH3 


25°S 


0.27 


z/4 


1080-1140 


0.29 


0.30 1.1 2.2 


H2S 


25°S 


0.27 


1^2 


1200-1250 


0.64 


91 330 700 



* Noise Equivalent Spectral Radiance averaged over the spectral interval, in 
units of nW cm~^ sr~^/cm~^. 
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Figure .1: (a): CIRS limb average spectrum of Titan from T55 flyby (black) and model 
fit (green), (b): data-model spectral residual. Note the I'j band of CaHg at 1158 cm"-'^ 
that is not included in our model (see text for details), (c) and (d) show the density and 
strengths of lines of PH3 and H2S respectively in the interval considered. 
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(a) T55 data and PH3 model 
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Figure .2: Calculation of upper limits for the abundances of H2S and PH3 in Titan's 
stratosphere, (a) & (c): residual of fit to the CIRS Titan spectrum in each region, after 
the modeling and removal known gas species (black line). Over-plotted is a calculation 
with exaggerated trial gas abundances of the undetected species to show their spectral 
signature (colored line), (b) & (d): the curve of growth of Ax^ over a wide range of trial 
abundances. The 1, 2 and 3-cr upper limits are indicated by the vertical dashed lines at 
Ax^ = 1, 4, 9 respectively. 
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